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We proposed a simple Ansatz for the vacuum wave functional (VWF) of SU(2) gauge theory in 
temporal gauge. In 2 + 1 dimensions, the Ansatz was shown to be a fairly good approximation to 
the true VWF of the theory. Relative probabilities of various test configurations in the vacuum can 
be computed in numerical simulations of lattice-regulated SU(2) gauge theory by the method pro- 
posed long ago by Greensite and Iwasaki. We report promising (albeit still preliminary) results of 
testing the proposed VWF in 3 + 1 dimensions on various sets of lattice gauge field configurations. 
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1. Introduction: An Ansatz for the vacuum wave functional of SU(2) Yang-Mills 
theory in temporal gauge 

The problem addressed in the present contribution is quite old, having been studied since at 
least 1979 |JJ]. It can be stated very simply: In the Hamiltonian formulation of SU(2) Yang-Mills 
(YM) theory in temporal gauge in D = d + 1 dimensions, one tries to find the ground-state (vacuum) 
wave functional (VWF) satisfying the YM Schrodinger equation plus the GauB law constraint: 

jH> [A] = J ^x|-i^^ + i^(x) 2 |^ [A] = E V Q [A], (1.1) 

^d k + ge^4) -J^VolA] = 0, (1.2) 

hoping to be able to extract from the VWF useful information on color confinement and relevance 
of various gauge-field configurations in its mechanism. As simple as it appears, the problem defies 
satisfactory solution, despite attempts using various techniques 1 (to appreciate the complications 
involved, see e.g. the horrifying expressions one deals with when looking for a solution perturba- 
tively [§). 

Some time ago, we proposed a simple Ansatz for an approximate VWF in temporal gauge [Q] : 



\Jd d xd d yFl){ 



1 




(1.3) 



In this expression, F?j{x) are components of the chromomagnetic field strength, S) 2 denotes the 
covariant laplacian in the adjoint representation, Ao is its lowest eigenvalue, and m is a parameter 
with the dimension of mass. 

In D = 2 + 1 spacetime dimensions, the above VWF is particularly simple. There exists 
a single magnetic field with 3 color components, B a {x) = F" 2 (x), a = 1,2,3, the square of the 
coupling constant g 2 has the dimension of mass, and the parameter m is to be chosen proportional 
to g 2 . We have shown that the proposed VWF Ansatz is a fairly good approximation to the true 
ground state of the theory. Analytical and numerical evidence in its favor is briefly summarized in 
Sect. | 



We conjecture that the expression ([L3J) could be a reasonable approximation to the true VWF 
also in D = 3 + 1 dimensions. However, to justify this proposition is much more complicated. 
Fortunately, there exists a method of measuring the (square of the) VWF in numerical simulations, 
proposed long ago by Greensite and Iwasaki [Qj, and utilized recently also for testing the wave 
functional in D = 2 + 1 dimensions [Q]. The method is outlined in Sect. [3] and the first preliminary 
results for test configurations in D = 3 + 1 dimensions are presented in Sect. |[ They are promising, 
but conclusions (Sect. ||) have still to be taken with some caution. 

2. Summary of results for (2+1) dimensions 

By construction, the VWF of Eq. (|0|) in (2 + 1) dimensions reproduces the well-known 

'Some approaches to the problem are reviewed in [0], references to original papers can also be found therein. 
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solution of the Schrodinger equation for the ground state of electrodynamics in the free-field limit 
(g — > 0). In addition, we have shown that it has a number of other attractive features: 

1 . It is a good approximation to the true vacuum also for strong fields constant in space and 
varying only in time [Q. 

2. For slowly vaiying chromomagnetic fields it reduces to the "dimensional-reduction" form: 

^o-exp -%jd 2 xB« l0 ^x)Bl lm {x) (2.1) 



that has been argued to be the correct form of the Yang-Mills VWF at large scales [|lj, 

3. It exhibits confinement (non-zero string tension) for the parameter m / 0, and the non-zero 
value of m seems preferred in a variational estimation of the vacuum energy [f|]. 

4. We proposed a recursion method of generating field configurations distributed according to 
the proposed VWF [Q] and compared a few quantities computed in an ensemble of recursion con- 
figurations with those obtained in ensembles of true SU(2) gauge-field configurations. The values 
of the mass gap the ghost propagator in Coulomb gauge, and the color-Coulomb potential [[7|] 
in recursion configurations reasonably agreed with true configurations. 

5. We computed the square of the true VWF for sets of simple test configurations (by the 
method described in some detail in Sect. |3|), the results were consistent with expectations based on 



the proposed VWF, Eq. (|L3j) 



The above results led us to conclude that the proposed VWF, in (2 + 1) dimensions, agrees 
with the true Yang-Mills vacuum wave functional for the bulk of the probability distribution, with 
possibly a small disagreement in the tail of the distribution. 

However, we live in D = 3 + 1 spacetime dimensions, and the success of our wave functional 



(1.3) in the real world is by far not guaranteed. If we postulate the same form to work also in 
D = 3 + 1 dimensions, there are some good and some bad news. The VWF again reduces to the 
QED VWF in the g — > limit (assuming m can be neglected as well), and also points 1 and 2 
above still remain valid [Q]. Unfortunately, the recursion method invented for generating gauge 
configurations in D = 2 + 1 does not work in D = 3 + 1. It profited from the possibility of easily 
going from potentials A* to field strengths B and back (in a variant of axial gauge), which is not 
possible in D = 3 + 1 because of Bianchi constraints. We are forced (at least for the moment) to 
test the proposed VWF only in numerical simulations on some sets of test configurations. 



3. Direct measurement of the VWF in numerical simulations 

There exists a method of direct measurement of the square of the VWF for some sets of gauge- 
field configurations in numerical simulations of the D = 4 Euclidean pure Yang-Mills theory. The 
method was invented by Greensite and Iwasaki [0], and is based on the observation that the square 
of the VWF for a certain gauge-field configuration U'(x) at fixed time (say t = 0) is expressed 
through a path integral: 

\V[U']\ 2 = ~J [DU]8(U )Ud [U(x,0) - U'(x)] e- s . (3.1) 

Now take a set of configurations ^ = {U {j \x),j = 1,2, . . . ,M}. The value of each |*F [U^] \ 
is proportional to Pj, the probability that a lattice configuration on the t = time slice is equal to 
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the 7-th configuration U^\x) £ . That probability can be computed numerically by a modified 
Monte Carlo simulation. All links at t ^ are updated in the usual way, e.g. by heat bath for 
the SU(2) Wilson action. At t = 0, one configuration from the set % is chosen randomly, and 
accepted/rejected by the usual Metropolis algorithm. Then 

pF[t/«]| 2 oc Mm (3.2) 



AW-*-" Ni 



tot 



where Nj is the number of times that - in a given simulation - the j-th configuration in the set is 
selected by Metropolis, and A^ot denotes the total number of updates of the t = plane. 

A limitation of the present method is that |*Po ] | falls °ff exponentially with the action 
of the configuration, so to achieve reasonable acceptance rate in Metropolis accept/reject steps all 
configurations in a test set have to be relatively close in action. (This is regulated by parameters K, 
a and y in sets of configurations below.) 



4. Preliminary results for test configurations in (3 + 1) dimensions 



In this pilot study, we address a simple question in lattice simulations: Can the data distinguish 
between the dimensional-reduction form of the squared VWF 2 



exp 



H d\B a k (x)Bl(x) 



(4.1) 



and that following from (1.3) 



exp 



ab 



d 3 xd 3 y B a k (x)\ 



B\{y) 



(4.2) 



What fits better to the true vacuum of the pure YM theory? 3 

To investigate this question, we computed, by the method described in Sect. ||, the squared- 
VWF values for three types of configurations: 

1. non-abelian constant configurations; 

2. abelian plane waves on lattices of variable size with maximum wave-length (A = L); 

3. abelian plane waves on a fixed-size lattice with variable wave-length. 

Simulations were performed for coupling constants j3 = 2.2,2.3,2.4, and 2.5 and lattices of sizes 
L = 12, 16,20,24, and 28. We present only a subset of data; more results with higher statistics and 
systematic error analysis will be presented elsewhere [gj. 



Non-abelian constant configurations: We used configurations 



NAC 



U k {j \x) = J\-(aU)y\ 2 + ia^o k 



with a 



(./) 



"if J- (43) 



2 On a lattice, the integral in the exponent is proportional to 5 sp = Y,p (l ~ ^Tr[t/f]), where the sum runs over all 
spatial plaquettes. 

3 The true VWF cannot just be of the dimensional-reduction form for all field configurations. This would lead to 
exact, not just approximate, Casimir scaling, and incorrect results at short distances/high frequencies | 
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for which 5 sp = K j. These are not useful to address the above question, because of the (approxi- 
mate) equality, valid for constant configurations: 




ab 



Bi(x 



1 _ I B b k (y)^-S sp . 



(4.4) 



(Sum over repeated indices is implied.) We used therefore these sets of configurations only to 
"calibrate" our computer code, i.e. for comparison with earlier data of Greensite and Iwasaki, which 
used modest-size lattices, 8 4 and 10 4 . New data, from 16 4 and 20 4 lattices, are fully compatible 
with the old ones. 



Abelian plane waves: The simplest choice is 



APW 



U[ j) (x) = \ll-[ d£> (x) ) 1 2 + ia [ n J > (x) a 3 ; U\ ]> (x) = U"> (x) = 1 2 



rU), 



(4.5) 



with 



a (j) (x) 



0C n + Ynj 



L 3 



■cos 



and X n = — . 

n 



For such configurations one expects 



Ni 



N u - 



(4.6) 



(4.7) 



where 



G>(Xn) 



from DR \¥ Q \ l , Eq. (4.1), 



from our V , Eq. (4.2), 



mr 



1-cos^^ . 



(4.8) 

At a fixed coupling j8, fixed lattice size L and a certain wavelength of X n , —log (Nj/N tot ) is 
a linear function of j, and one can get (o{X n ) from the best fit of the form (fOJ). In principle the 
obtained slope ft) could depend on the choice of constants a„ and y n that parametrize the set of 
configurations. It turns out that slopes coming from quite different pairs (a„,y„) are compatible 
within errors. This is illustrated in Figure |] One should however keep in mind that the ambiguity 
of the choice of (cc n ,y n ) is a source of systematic error that has not yet been accounted for in the 
analysis of data presented below. 

Figure ^ represents the results for abelian plane waves with maximum wavelength X = L 
for a number of lattices sizes L at a fixed gauge coupling, j8 = 2.4. The solid (red) line repre- 
sents a fit of the form x k 2 (L) / ^k 2 (L) +m 2 ^ motivated by our proposed VWF The data are 
described satisfactorily by this fit. However, even the fit of the form (a + bk 2 {L)^ (motivated by 
the dimensional-reduction VWF) does represent the data quite reasonably (the dotted blue line), 
though the % 2 of the latter fit is considerably worse than of the former one. 

The parameter m of the former fit, if it corresponds to a physical mass in the continuum limit, 
should scale properly as a function of the coupling j3. The preliminary data look promising in this 
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p=2.5, 28 4 lattice 



1 — H — 1 par. set 1 




- slope=0.175(1) 




■ par. set 2 




- slope=0.181(4) 






* 




p 








w 








/ 











20 40 60 80 100 120 140 
0.5(a + H) 

Figure 1: The slope of the linear dependence of — log (Nj /Mot) vs. \ {oL + yf) depends slightly on the chosen 
parameter set (a, 7). The displayed data are for abelian plane waves with A = L at /3 = 2.5 and L = 28. 



respect. If we divide m by the value of another mass quantity, e.g. the square root of the string 
tension, at the same coupling, the ratio should not depend on the coupling. In Figure |3| we display 
such ratio at four values of gauge coupling in the scaling window, and the result is compatible with 
a constant, m/y/o ~ 2.36. One should, however, take this result with some caution: the errors 
shown may be underestimated, a substantial change of m can be compensated by a change of the 
fit parameter c without deteriorating the agreement with data. This is illustrated in Figure [| A 
fit to data with one parameter set (a, 7) at j8 = 2.2 gave m = 1.45(40), a different parameter set 
(a, 7) leads to m = 1 .07(8), and the fit to combined data (that was displayed in Figure ^ resulted in 
m = 1.17(14). Still, the two fits with m differing by ~30% are hardly distinguishable in Figure |J 
To discriminate between the two forms in Eq. ( |4.8[ ) one needs to use sets of abelian plane 
waves with variable wavelengths X n on a fixed size lattice. An illustrative result is displayed in 
Figure ||. Here the fit motivated by our proposed form of the VWF is clearly preferred. However, 
the mass parameter m resulting from the best fit to these preliminary data (m w 1.29) is somewhat 
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Figure 3: The ratio of the parameter m (extracted from fits to data for abelian plane waves with maximum 
wave-length) to the square root of the string tension a vs. the coupling constant /3. 

different from the value obtained from abelian plane waves with maximum wavelength (m 0.90). 
More extensive simulations are needed to resolve this discrepancy. 

5. Conclusions 



The preliminary investigation of the SU(2) Yang-Mills vacuum wave functional in D = 3 + 1 
dimensions has provided promising outcome: 

1. The results of this pilot study are consistent with the old ones by Greensite and Iwasaki [9], 
but extend their investigation to larger lattices, and direct tests of the approximate form of the 



vacuum wave functional ( \l.3\ ). 

2. For the three sets of simple configurations (non-abelian constant configurations, abelian 
plane waves with X = L, abelian plane waves with variable wavelengths) the computed values of 



|*Fo| 2 are consistent with the form (1.3), but the evidence cannot be considered conclusive. 
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Figure 4: co vs. k(L) for abelian plane waves with maximum wave-length at j3 = 2.2 for a number of lattice 
sizes L. Two fits to data sets corresponding to different choices of parameters (a, 7) give considerably 
different m values. 
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Figure 5: (0 vs. k(X) for abelian plane waves with variable wave-length at j3 = 2.3 on 24 4 lattice. 

3. The parameter m of the VWF, determined from fits to probabilities of plane- wave configu- 
rations with X = L, scales as a physical quantity with dimension of mass. 

4. The numerical data for abelian plane waves deviate from expectations based on the dimen- 
sional-reduction form of the VWF. 

For the future, more extensive simulations, higher statistics and systematic error analysis for 
simple test configurations are called for. The most important (and difficult) task is to find an effi- 
cient method in (3 + 1 ) dimensions for generating gauge-field configurations distributed according 
to the VWF dOl ). 
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